On the enstrophy dissipation in two-dimensional turbulence 
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Insight into the problem of two-dimensional turbulence can be obtained by an analogy with a heat 
conduction network. It allows the identification of an entropy function associated to the enstrophy 
dissipation and that fluctuates around a positive (mean) value. While the corresponding enstrophy 
network is highly nonlocal, the direction of the enstrophy current follows from the Second Law of 
Thermodynamics. An essential parameter is the ratio = ~fk/(vk 2 ) of the intensity of driving 
7fc > as a function of wavenumber k, to the dissipation strength vk 2 , where v is the viscosity. 
The enstrophy current flows from higher to lower values of Tfe , similar to a heat current from higher 
to lower temperature. Our probabilistic analysis of the enstrophy dissipation and the analogy with 
heat conduction thus complements and visualizes the more traditional spectral arguments for the 
direct enstrophy cascade. We also show a fluctuation symmetry in the distribution of the total 
entropy production which relates the probabilities of direct and inverse enstrophy cascades. 
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I. INTRODUCTION 

Three-dimensional turbulence displays an inertial 
range, in which energy is transferred from the spatial 
scales at which it is introduced into the system down 
to small scales, where it is finally dissipated by viscous 
forces. The standard picture of turbulence in two dimen- 
sions is qualitatively different. Following the pioneerin 
works of Kraichnan p], 0, Leith and Batchelor 
(KLB) on the two-dimensional Navier-Stokes equation 
for the fluid velocity field, one expects an inverse energy 
cascade from the forcing scales to large scales and simul- 
taneously a direct enstrophy cascade from the forcing 
scales to small scales. The enstrophy is the variance of 
the vorticity, namely the ensemble average of the squared 
curl of the velocity. 

Two-dimensional turbulence has been a very active 
area of theoretical, numerical and experimental inves- 
tigation 0, , not only as an easier test case but also 
relevant to certain real quasi-two-dimensional situations. 
Examples include oceanic currents and atmospheric and 
geophysical flows 0, but two-dimensional flow is also 
realized in laboratory situations 0, @ . However, the pic- 
ture of two-dimensional turbulence remains not fully un- 
derstood and in fact, there are some limitations to the 
above classical KLB scenario. For example, the stan- 
dard enstrophy cascade disappears when considering a 
bounded domain where only a monoscale forcing is ap- 
plied Moreover, the mechanism of the direct enstro- 
phy cascade and the determining factor for the direction 
of the enstrophy current has not been fully understood 
as a consequence of a more general principle. There have 
been recent clarifications, g oing into the details of the 
physical mechanism, e.g. [j| [H], but it seems interest- 
ing and natural to connect the situation also with better 
understood scenario's and to be able to see the enstro- 
phy dissipation as the result of a more generally valid 
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principle. 

In the present paper, we address the issue of the en- 
strophy current and of its direction. A very close anal- 
ogy with a two-dimensional heat conduction problem 
provides new ingredients to understand the enstrophy 
cascade in its full qualitative behavior. It turns out, 
as will be shown later, that the stochastically driven 
Navier-Stokes equation for the vorticity can be mapped 
to a problem of heat conduction: at each wavenum- 
ber k a thermal reservoir is attached with temperature 
Tk = jk/ivk 2 ) where 7& is the forcing strength and v is 
the viscosity. From the Second Law of Thermodynam- 
ics, that will be derived in its detailed version, follows 
that enstrophy is dissipated as heat flows: from higher 
to lower temperature, or here, when 7^. is peaked around 
some small mode fc, from small to large wavenumbers. 
In other words, the origin and the direction of the en- 
strophy flux is simply and directly a consequence of the 
Second Law applied to the enstrophy. We also go beyond 
the study of the average enstrophy current and discuss a 
symmetry in its fluctuations. That estimates the prob- 
ability of going backwards, i.e., the probability of an in- 
verse enstrophy cascade. At the same time, we obtain for 
the first time a steady state fluctuation theorem in the 
context of turbulence. 

In the next section, we start by reminding the reader 
of the standard picture of two-dimensional turbulence. 
In Section IIIII comes the analogy with heat conduction. 
From it follows the final analysis of the enstrophy dis- 
sipation in Section IIVI The main general consequences 
and conclusions are taken in Section Ivl 

The paper will describe the arguments and analogues 
in a formal way, avoiding however a fully rigorous math- 
ematical analysis. The main goal is indeed to point out a 
useful picture and analogy which is sufficiently powerful 
to specify the enstrophy cascade. To add the mathemati- 
cal details and hypotheses is not believed to be extremely 
difficult but only few remarks are added to guide the 
mathematically inclined. 
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II. NAVIER-STOKES EQUATION 

The Navier-Stokes (NS) equation [H E2 for the ve- 
locity field u(t, r) is 

Oil -> -* -> 

— + (u-V)u = vAu-Vp + f (III) 

where p is the pressure, / is the external force, and v is 
the viscosity. That is supplemented by the incompress- 
ibility condition 

V-u = 0, 

and in our case by periodic boundary conditions for a 
finite spatial region V. Similar equations arise for the 
vorticity lj = V x u, by taking the curl of 

^ + (u ■ V)lo - (Q ■ V)u = i/A5 + g (II.2) 

with g — V x f . The energy of the system is given by 
the total kinetic energy 

£ = Ay 

while the enstrophy is defined as 

Jv 2 

Its role will become clearer later on. 

Consider Cartesian coordinates r = (x±, X2, X3) = 
(x,y,z); the two-dimensional case is conveniently rep- 
resented by setting the third component of the velocity 
equal to zero: u— (ui,u 2 ,0). Therefore Q = (0,0, lu 3 ) is 
better represented by a pseudo-scalar u = 0J3. Eq. illl.2|) 
thus acts on a single component and since u and Q are 
now perpendicular to each other, it is further simplified 
by the vanishing of the term (uj ■ V)w, 

du> - 

— + (u ■ V)uj = vAuj + g (II.3) 

where g = dfi/dy — df2/dx. The pressure has dis- 
appeared but equation (|II.3|I is still nonlocal because 
u = K,uj for some Biot-Savart kernel IC. 

We take our system bounded in a rectangular domain, 
where it is useful to consider the Fourier transform 

These modes satisfy uj k = uj- k which will always be un- 
derstood. 

Upon Fourier-transforming 1|II.3|) we thus get, for k £ 
Z 2 \{0}, 

^~F k (Lu) = -vk 2 Lu k +g k (II.4) 



where 

F k (uj) = QjjLUjUjf (II. 5a) 

j+e=k 

with coefficients given by 

These <F,,£ = if and only if either j \\ £ or \j\ = \£\. 
Alternatively, 

Fk(u) = 4> k ,i uiujk-e (II. 5b) 

with 

, _ h£ 2 - hk 2 (I 1 \ 

^ e = — ^ — [w Jk~W ) ^ e 

In that notation, 

(j>k,e + <Pk,k-e = 

represents the so called triad relation of 13] . Another 
alternative, which will be useful later, is 

F k(") = -^1^ \ k _ e \2 Wk-t • (II.5c) 

Finally we must specify the forcing g k . A transla- 
tionally invariant and stationary turbulent state can be 
achieved by imposing a force that is homogeneous in 
space and time. A Gaussian random field with zero mean 
is the simplest example: in that case the force / in 1II.1|) 
is a Gaussian noise that is white in time and colored in 
space, completely determined by its covariance 

(f i (s,r)f j (t,r'))=C ij (r-r')S(t-s), 

where diCij = (incomprcssibility). Equation 1II.4|) 
then turns into the stochastically driven NS 

duj k {t) = -vk 2 Lo k dt + F k {uj)dt + v / 2r/ k ~dW k (t) (11.6) 

in which dW k — dW~ k represents a standard Wiener 
process |l4j . That driving pumps vorticity into the sys- 
tem at wavenumber k with intensity ^ k > while the 
viscosity v > enters in the first term on the right-hand 
side of (|II.6(1 to dissipate the vorticity. Equation (|II.6(I is 
the starting point of our analysis. 

A. Mathematical assumptions 

The point of departure (|II.6|I is a stochastic differ- 
ential equation to be understood in the Ito-sense [l4| . 
Solutions are Markov processes but note that they are 
infinite-dimensional. In general the resulting diffusion is 
not elliptic because some 7^ can be made zero. That 
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brings us to the problem of understanding the assump- 
tions on the strengths 7^ and on the viscosity v so that 
there is a unique invariant probability measure fi. A lot 
of mathematical work has been devoted to that problem 
in recent years. For example |l5j . if v is sufficiently large 



and 



as a function of J2"/k> tnen A 4 ^ s unique. Also \i is 



unique when there is k > so that 
a and for every |fc| > k. Even |l7l IT 



|fc| a for some 
)], when 7fe ^ 



for every \k\ < N where N is some number that depends 



epe i 

on v and on X) 7fe: then /1 is unique. We refer to p0| for 
even stronger and more recent results. 

In what follows we simply assume, with no further 
ado, that fi is unique and has smooth local densities. 
Another assumption in the technical manipulations is to 
start from a finite dimensional analysis. In other words, 
we choose a finite but arbitrary N and consider equation 
jlL6)l only for k 2 < N with j k f there. That cut-off 
will take care of convergence problems in what follows 
and it allows us to speak of fi(u>) as the density of /x with 
respect to the flat measure du. 



B. Euler equation 

The vorticity and the corresponding enstrophy play an 
important role in two-dimensional turbulence because of 
the appearance of an extra conservation law. The Euler 
equation corresponding to (|II.6(1 is 



du> k {t) = F k (uj) dt . 
cf. (|II.5a|l l|II.5c(l . It is easy to see that 

k 

so that enstrophy is conserved 

dQ v l d\u; k \ 2 



(II.7) 



2 dt 



As a consequence, the enstrophy in (|II.6J| is changed by 
the injection (at rate 7^) and the dissipation (with in- 
tensity v) but is transported without dissipation over 
the various modes via the nonlinear and highly nonlo- 
cal terms in F k . That invites the definition of various 
enstrophy currents. 



C. Currents 

The net enstrophy current J k that leaves the system at 
wavenumber k is obtained from investigating the sources 
and sinks to the enstrophy. The total enstrophy dissipa- 
tion over the time interval [— r, r] is computed from 



ft(t) = ~$> fc | 2 (i) 



\oj k \ 2 {t)dt- 



2 7fc 5i/ W k {t)odW k (t) 



where the last integral is in the Stratonovich sense [l4| . 
and 3? stands for real part. The expression thus evaluates 
the change of enstrophy during a time interval [— r, r] 
for a history (io k (t)) and for a realization of the noise 
(dW k (t)). It is therefore natural to put 



TOUt 
J k 



vk 2 J \uj k \ 2 dt 
J k n = ^271 » / u k o dW k 



as the current going "out" , respectively "in" the system 
at mode k with respect to the external enstrophy reser- 
voir. The difference, 

7 rout rill 

Jk = J k - J k 

is the net enstrophy current that leaves the system (en- 
ters the environment) at mode k. 

On the other hand the local conservation law reads 



K| 2 (r) M 2 (-r) 



= J \uk\ 2 (t)dt 

+5R J uJ k F k {io)dt 
+ ^27^ / uj k odW k (t) 



which defines 

Jlk = 



2k 



dt 



hk 2 - h£ 2 
\£~k\ 2 



■ w-k^k-e^e 



the net current from mode I to mode k [here we used 
(|11.5c[) ] . Note the asymmetry J k g — —Ji k . 

As said before, the redistribution of enstrophy due to 
interactions between different modes globally does not 
change the total amount of enstrophy in the system: 
E k J2 e Jki = 0, see(|lL7|>. 

One of the main problems for the cascade picture is to 
understand the direction of the flow of the J k g. That is 
basically determined by the stationary (J k ). At the end 
of the paper we also discuss its fluctuations. 



D. Spectral distribution 

Heuristically, the reason why the enstrophy flows to- 
wards small scales (large wavenumber k) is because at 
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these small scales the dissipative term v/\uj in l|II.3|) dom- 
inates over the advection term (u ■ V)u. A more refined 
argument, started by Kraichnan derives the two- 
dimensional cascade picture (the so called direct cascade 
for the enstrophy and the inverse cascade for the energy) 
by investigating the energy spectra. The Fourier spec- 
trum of energy embodies the KLB picture by showing a 
power-law regime for each of the two cascades. Since the 
enstrophy spectrum is simply related to the energy one, 
from the inspection of the energy spectrum one can argue 
where energy and enstrophy are transferred or dissipated. 
In a way the cascade of enstrophy to small scales is the 
two-dimensional analogue of the energy cascade in three 
dimensions. We skip here the details of that Kolmogorov- 
Kraichnan theory as t hey have been excellently reviewed 
by many. We refer to |2l], |23| ■ 

The purpose of the present paper is to try an alter- 
native to that analysis by opening an analogy with heat 
conduction. It is interesting that in this way a natural en- 
strophy dissipation function appears, the thermodynamic 
entropy production, as will be explained in Section llVI 



FORMAL ANALOGY WITH HEAT 
DISSIPATION 



III. 



Remember our starting equation i|II.6|) . Let us first 
forget about the coupling between the various modes so 
that the system is reduced to the stochastic dynamics 



duj k {t) = ~vk 2 uj k {t)dt + ^j2^dW k {t) 



(III.l) 



describing an ensemble of uncoupled oscillators labeled 
by the wavenumbcr k. While in the original NS equa- 
tion the viscosity represents an irreversible loss, here it 
balances reversibly with the stochastic forcing. The dy- 
namics Ijlll.lt has a reversible equilibrium measure 



M°(^)=n 



\ui k \ /2T k 



■ dw k doj k 



(III.2) 



that we will use as a reference. 

The parameter T k = jk/ {vk 2 ) can be viewed as a kind 
of "temperature" of the reservoir attached to wavenum- 
ber k; it is of course no physical temperature. Thus, 
our approach is different from previous attempts to use 
a thermodynamical formalism in turbulence, identifying 
variables like io\ with a temperature (see (23| and refer- 
ence therein). 

The reversibility of the dynamics l|III.l|) is taken with 
the usual kinematical time-reversal that reverses the sign 
of the velocity field: the dynamical time-reversal of a 
history 

£=(w(*),*e[-T,r]) (III.3) 
in a given time interval [— r, r] is 

e£ = (-w(-t),te[-T,T]) (III.4) 



When we add the F k (ui) to (|III.1(I to obtain HII.6fl the 
oscillators become coupled, in fact in a nonlinear and 
nonlocal way. That coupling does however preserve the 
enstrophy very much like a Hamiltonian coupling that 
conserves the energy. The picture that thus emerges is 
formally equivalent to a heat conduction network where 
the vertices of the network are represented by the modes 
k. 

The Euler equation represents the conservative part of 
the time-evolution. That is changed by the addition of 
the Langevin forces that represent "thermal" reservoirs 
at each of the k; thus obtaining our equation l|II.6|) . Ob- 
serve that the "friction" depends on the "location" k of 
the oscillator. Standard thermodynamics then teaches 
us that there will be a "heat current" from higher to 
lower temperature. That "heat current" is in our present 
set-up played by the enstrophy current. Hence, if the 
driving makes T k a decreasing function of |fc|, e.g. by 
having 7^ ~ k~ a for some a > 0, then, the enstrophy 
should be transported from small |fc| towards larger \k\. 
That picture will be detailed in the following sections. 

The forward generator C + corresponding to the 
Markov diffusion 1II.6|) can be split into a "conservative" 
and a "dissipative" part, C + = + C\, with 



c-tp 

and 



E 

k 



F k {to\ 



dp_ 

duj k 



(III.5a) 



duj k 



y 



dXk 
du>k 
(III.5b) 



where we made use of the shorthand 



X = e -/3 fc K| 2 /2 



d_ 

dujk 



(e 



/3fcK| 2 /2 



p) 



with 0k = vk 2 /-f k . 

For the stationary measure fi we have = and in 
particular 



{F k {uj)uj k } - vk 2 {\uj k \ 2 ) +7 fc 

From now on we use that notation < 
tionary average according to [i. Equation l|III.6|) gives, 
for every time interval [— r, r], 



. (III.6) 
to denote a sta- 
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(J fe ) =vk 2 {{\u k \ 2 ) -T fc ) 



(III.7) 



That equation is the detailed enstrophy balance equation 
in stationarity; summing over k gives the somewhat more 
familiar 



k 



but at the same time and as a new interpretation of 
HIII. 70 we recognize how the net current into the enstro- 
phy reservoir at mode k is like a heat current into a ther- 
mal reservoir as determined by the difference between, 
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what now plays the role of a local kinetic temperature, 
(|wfe| 2 ) and the reservoir temperature T k . 

IV. ENSTROPHY DISSIPATION 

Continuing with the analogy above a quantity is now 
brought to the forefront which we call the entropy current 
S. Since the net enstrophy current leaving the system at 
each mode k is J% — J k ut — J™ and the corresponding 
"effective temperature" is T k we put 

S^-U (IV.l) 

h k 

as variable entropy current. It is a function of the history 
(|III.3|I over [— r, r]. The entropy current S is the entropy 
production in the environment associated to the enstro- 
phy dissipation; it is the usual sum over all dissipative 
currents divided by the respective temperatures. In the 
stationary state, the average (S) is the total change of 
the entropy in the universe over the time-interval [— r, r]. 

We will now show what is suggested thermodynami- 
cally by the previous analogy: S should measure the ir- 
reversibility and {S) > as a consequence of the Second 
Law of Thermodynamics. 

Remember that fi is the stationary measure of the NS 
dynamics (|II.6(1 : we denote by 7r/x its time-reversal. In a 
given time interval [— r, r] each history (|III.3() is realized 
in the system with a probability that comes from the 



path-space measure PZ(d£), i.e., the stationary Markov 
diffusion process associated to the stationary measure \x 
and the stochastic dynamics 6|) . 

We compute the logarithmic density (see also <|III.3() - 

X^llh (IV - 2) 

as a measure of irreversibility. It gives the ratio between 
the probability of a history £ and the probability of the 
time-reversed history 0£. We show that R coincides with 
S up to a temporal boundary term. Moreover, taking 
stationary averages (R) = (S). Since by construction, 
(R) > it also follows that (S) > 0. 

To compute R it is useful to compare the path-space 
measure with the reference path-space measure of the 
uncoupled case (|III.1|> . denoted by P®f (that one is sta- 
tionary and reversible). Thus first we compute the action 

pr 

V 

and similarly A^^ o 0, to finally estimate (|IV.2J| as the 
source of time-reversal breaking 

R = - Ar M o 6 (IV.3) 

The comparison of the two measures P and P° is made 
by means of the Girsanov formula 01 > obtaining 



V — 



V k 2 ^.{uj k F k {uj)) - -\F k {u)\ 2 dt + Vt 



F k (uj)du! k 



+ ln M (w(-r)) -hxfJP(uj(-T)) 



Substituting 0£ gives 



-vk 2 $t(u k F k (u;))-^\F k (u)\ 2 



F k (uj)duj k 



oe|+ln/i(w(r))-ln/x°(a;(T)). 



Here Ito stochastic integrals are performed and one should remember that these are themselves not time-reversal 



symmetric 14] . As an example for computing l|IV.3|) we see that 



lim 

At^O 



^2 F k (u {tj-i )) [uk (tj )-u k (tj-i)\ 



lim J? 

At^O 



F k (uj)dui k 



F k (u)duj k 



o6 



Fk (-wfe)) [-w fc (t j-i ) + uj k (t j ) 



lim 3? 

Ai^O 



W k (t 3 ) — UJk(tj-l) 



T dF k (u) 

-r d^k 



dt 



= 



because dF k /du k = 0, see (|II.5b|) . 



As a consequence, (|IV.3|) becomes 

R = S + In h{w{-t)) - In h{u{t)) 



with 

(IV.4) 5 = V — 
^T k 



%l[Lu k F k (uj)}dt 
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From 8(1 that expression coincides exactly with (|IV.1(I , 
as promised. 

When instead of the stationary p we had taken some 
initial density evolving as pt,t G [— r, r], the analysis 
above would be essentially unchanged. In that case the 
source of irreversibility is 



R = ^ 7-7- Jk + hxp- T (u- T ) - In p T (v T ) 

-L k 



(IV.5) 



where the only difference with (|1V.4|I resides in the last 
two terms, the temporal boundary 

[-lnp T (uj T )] - [-lnp_ r (w_ r )] 



A. Mean entropy production 



From the definition ljiV.2|) it directly follows 

(e"«) = 1 

(it is essentially the normalization condition of the path- 
space measure P^Q). Hence, by a convexity inequality, 
the stationary enstrophy dissipation (S) = {R) > 0. We 
can however be more explicit concerning that point by 
deriving an expression for (S) which is explicitly non- 
negative. In fact, we will show that 

(S) = E7 fe ((exphy fe H]-^exp[^H]) 2 } (IV.6) 



where 



Vfc(w) = K| 2 /(2T fc )+ln M (w) 



From (|IV.6fl . (S) > strictly as we can only have 

d 



dtok 



/3 k \u> k \ 2 /2 



p{uj) 



= 



for all k when p = po of ljIII.2(l . 

Here comes the proof of (|1V.6|) . Denote by E p _ T the 
expectation in the process P^_ x started from p_ T . We 
assume that at time r the evolved measure is described 
by a density p T . We have (|IV.5(I . in expectation, 

E P- T [R] = J2 I J fe] + S (Pr) - S (P-r) (IV.7) 

k 

where S(p) = — J duj p{oj) In p(u>) is the Shannon entropy 
of the density p. Another formulation is 



,-r[R] 



R(t) dt 



with, similar to (IIII.7I) 



R(t)=vY. k2 ^j du>\cJ k \ 2 p t (u)-j-] + j t S(p t ) 



(IV.8) 



The previous considerations thus identify the mean dissi- 
pation rate at time t (in the transient regime) with R(t). 

To see the relation with ljIV.6|) we start by evaluating 
the time-derivative of the Shannon entropy: 



dS_ 
~dt 



dm — In p 
dt ' 



dw(C + p)bxp (IV.9) 



Using the invariance of the Shannon entropy under the 
conservative (Euler) part of l|111.5a(l . we get 



duj k 



d ^(p) = - [ duj(C+p)\np = ^~i k [ duX k d] '' ! 

J k J 



Minus the second term reads 



l/^^k 2 J du uj k X k — v £ k 2 J du> u: k 

k k " 

= v^2f3 k k 2 J duj p[\uj k \ 2 - 



(IV.10) 



Pk^kP 



(IV. 11) 



Subst ituting IjlVTll) into (|IV.10jl and then (|IV.10jl int o 
(|IV.8(1 . we immediately obtain the desired identity (|IV.6|I . 



B. Enstrophy network 

The situation can now be summarized as follows: lo- 
cally, in the stationary measure, we have 



and globally 



£<-4> = o 



For the enstrophy dissipation 



(S) =Y J Pk(J k ) >o 



(IV.12) 



(rv.13) 



We have here formally the same situation as for a heat 
conduction network as considered e.g. in [24j . The rela- 
tions (|IV.12|I and (|IV.13|) do not of course uniquely deter- 
mine the mean enstrophy currents but their direction or 
sign is thermodynamically determined by analogy with 
heat conduction. 

Let us first consider the typical case where the 
strengths 7^ are non-zero only for a neighborhood of 
k = 0, say j k = 1 when |fc| < r\ and outside that 
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large wavelength regime, jf. I 0, \k\ > rj. In terms of 
heat conduction it would mean that the temperatures 
Tk = \/{vk 2 ) are decreasing outward in the disk for 
| As | < 77 and fall to Tk = outside (|A:| > rf). Clearly 
then, there will be a heat current toward increasing |fe| or, 
here, an enstrophy current towards smaller wavelengths. 
In other words, the enstrophy current is a kind of non- 
local heat current the direction of which is determined 
by the Second Law. Because of the nonlocality of the 
term Fk{uj) the current will not stop at the boundary 
of the disk but will be more and more suppressed when 
regarding Jki for k inside and £ outside the disk. For re- 
ally large i there is no longer a visible local heat current. 
That seems compatible with the observations [8j that the 
enstrophy cascade remains pretty localized around the 
forcing window. 

In general however, when all 7^ > are active, we 
have a "temperature" profile jk/i^k 2 ) that can of course 
be complicated. If the jk only depend on |As| we have 
in essence a one-dimensional heat conduction problem 
(along the radial direction). 



which is sometimes easier to check numerically and ex- 
perimentally. Roughly speaking, that last relation tells 
us that the probability of observing the inverse cascade 
for the enstrophy is exponentially smaller than the prob- 
ability of observing the direct cascade. 

V. CONCLUSIONS 



The main conclusion is derived from the analogy with 
a heat conduction network. Above and beyond all de- 
tailed physical mechanisms that give rise to the direct en- 
strophy cascade in two-dimensional turbulence stands the 
Second Law of Thermodynamics for the entropy l|IV.l|) 
which gives a direction to the enstrophy flow. The rel- 
evant parameter is the ratio ^k/ivk 2 ) which plays the 
role of an effective temperature of an enstrophy reservoir 
to which each mode k is coupled. If the forcing is re- 
stricted to a finite window, then the temperature outside 
is effectively equal to zero. The conservative part in the 
enstrophy conduction is nonlocal but does not contribute 
to the dissipation. 



Fluctuations 



Looking back at ijIV.7(l and l|IV.8|) , we found the mean 
entropy as the change of entropy in the environment S 
plus the change of (Shannon) entropy due to the stochas- 
tic dynamics in the system. Its stationary mean (R) > 
is strictly positive. We will now look at its fluctuations. 
More precisely, we consider the R of (|IV.4|) and ask for 
its probability distribution. Since by construction (|IV.2|1 . 



we have that 



J f(S0dP^ = J f(0e- R{i) dP^) (IV.14) 

is exactly valid for all observations / and for all times 
r. The relation ljIV.14(l is called a fluctuation symmetry, 
see e.g. j2^, because it generates the so called fluctuation 
theorem for the entropy production as first formulated 
in j2(| I27I l28l ] . Remember that R equals the S up to a 
temporal boundary term, see l|IV.4|) and l|IV.5|) . 

One of the consequences of the fluctuation symmetry 
(I1V .14| is that 



Prob[i? < 0] 
Prob[i? > 0] 



|i?>0> 



(IV. 15) 



We have identified a general entropy function l|IV.l|) 
and (|1V.4|I . also in the transient regime, see (|1V.7() and 
(|IV.8|I . We have shown that the stationary entropy pro- 
duction is strictly positive. It provides the general mech- 
anism driving the direct enstrophy cascade. The fluctua- 
tions in the entropy satisfy the symmetry l|IV.14(l which 
gives an estimate (|IV.15ll of the relative probabilities of 
direct versus inverse cascades. 

An important open question remains however. The 
above analogy is silent about the inverse energy cascade. 
We have not found a heat conduction analogue which 
would reveal the inverse cascade for the energy dissipa- 
tion in two-dimensional turbulence. Of course, as energy 
and enstrophy are entangled and spectrally related, the 
direct enstrophy cascade has direct consequences in the 
form of the inverse energy cascade. That point follows 
from the standard treatments, see also Section Hi Dl as 
in but has not been clarified in the present paper. In 
fact, a naive extension of the present formalism but for 
the energy would find the inverse cascade quite surpris- 
ing as it seems to reduce entropy. We have not inves- 
tigated whether the combination of dissipative currents, 
enstrophy and energy, would still lead to a total posi- 
tive entropy production, as expected thermodynamically. 
Clarifying that entropy balance remains one of the most 
intriguing problems of two-dimensional turbulence. 
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